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Abstract 

John Tromp introduced the so-called ’binary lambda calculus’ as a way to encode lambda 
terms in terms of binary words. Later, Grygiel and Lescanne conjectured that the number of 
binary lambda terms with m free indices and of size n (encoded as binary words of length n) 
Is o for r ~ 1.963448 .... We generalize the proposed notion of size and show that 

for several classes of lambda terms, including binary lambda terms with m free indices, the 
number of terms of size n is 0 p~"'j with some class dependent constant p, which in 

particular disproves the above mentioned conjecture. A way to obtain lower and upper bounds 
for the constant near the leading term is presented and numerical results for a few previously 
Introduced classes of lambda terms are given. 


1 Introduction 

The objects of our interest are lambda terms which are a basic object of lambda calculus. A 
lambda term is a formal expression which is described by the grammar M ::= x \ Xx.M \ {M N) 
where a; is a variable, the operation (M N) is called application, and using the quantifier A is called 
abstraction. In a term of the form Xx.M each occurrence of a: in M is called a bound variable. We 
say that a variable x is free in a term M if it is not in the scope of any abstraction. A term with 
no free variables is called closed, otherwise open. Two terms are considered equivalent if they are 
identical up to renaming of the variables, i.e., more formally speaking, they can be transformed 
into each other by a-conversion. 

In this paper we are interested in counting lambda terms whose size corresponds to their De 
Bruijn representation {i.e. nameless expressions in the sense of 0 ). 

Definition 1. A De Bruijn representation is a word described by the following specification: 

M ::= n| AM|MA^ 

where n is a positive integer, called a De Bruijn index. Each occurrence of a De Bruijn index is 
called a variable and each X an abstraction. A variable n of a De Bruijn representation w is bound 
if the prefix of w which has this variable as its last symbol contains at least n times the symbol X, 
otherwise it is free. The abstraction which binds a variable n is the nth X before the variable when 
parsing the De Bruijn representation from that variable n backwards to the first symbol. 

For the purpose of the analysis we will use the notation consistent with the one used in [T]. 
This means that the variable n will be represented as a sequence of n symbols, namely as a string 
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of n — 1 so-called ’successors’ S and a so-called ’zero’ 0 at the end. Obviously, there is a one to one 
correspondence between equivalence classes of lambda terms (as described in the first paragraph) 
and De Bruijn representations. For instance, the De Bruijn representation of the lambda-term 
Xx.Xy.xy (which is e.g. equivalent to Xa.Xb.ab or Xy.Xx.yx) is AA21; using the notation with 
successors this becomes AA((S'0)0). 

In this paper we are interested in counting lambda terms of given size where we use a general 
notion of size which covers several previously studied models from the literature. We count the 
building blocks of lambda terms, zeros, successors, abstractions and applications, with size a, b, c 
and d, respectively. Formally, 

|0|=a, \Sn\ = \n\+b, \XM\ = \M\+c, \MN\ = \M\ + \N\ + d. 

Thus we have for the example given above |AA((5'0)0)| = 2a -I- 6 -I- 2c -I- d. Assigning sizes for 
the symbols like above covers several previously introduced notions of size: 

• so called ’natural counting’ (introduced in [T]) where a = b = c = d=l, 

• so called ’less natural counting’ (introduced in [T]) where a = 0, & = c = 1, d = 2. 

• binary lambda calculus (introduced in |^) where 6=l,a = c= d = 2, 

Assumption 1. Throughout the paper we will make the following assumptions about the constants 
a, b, c, d: 

1. a,b,c^d are nonnegative integers, 

2. a -f d ^ 1, 

3. b,c> 1, 

4- gcd(&, c, a -I- d) = 1. 

If the zeros and the applications both had size 0 (i.e. a + d = 0), then we would have infinitely 
many terms of the given size, because one can insert arbitrary many applications and zeros into 
a term without increasing its size. If the successors or the abstractions had size 0 (i.e. b or c 
equals to 0), then we would again have infinitely many terms of given size, because one can insert 
arbitrarily long strings of successors or abstractions into a term without increasing its size. The 
last assumption is more technical in its nature. It ensures that the generating function associated 
with the sequence of the number of lambda-terms will have exactly one singularity on the circle of 
convergence. 

Notations. We introduce some notations which will be frequently used throughout the paper: 
If p is a polynomial, then RootOf {p} will denote the smallest positive root of p. Moreover, we will 
write [z"]/(z) for the nth coefficient of the power series expansion of f(z) at z = 0 and f(z) -< g{z) 
(or f(z) < g(z)) to denote that [z”] f{z) < [z”] g(z) (or [z"] /(z) < [z”] ^(z)) for all integers n. 

Plan of the paper. The primary aim of this paper is the asymptotic enumeration of closed 
lambda terms of given size with the size tending to infinity. In the next section we define several 
classes of lambda terms as well as the generating function associated with them, present our main 
results and prove several auxiliary results which will be important in the sequel. We derive the 
asymptotic equivalent of the number of closed terms of given size up to a constant factor. This 
is established by construction of upper and lower bounds for the coefficients of the generating 
functions. These constructions are done in Sections and To get fairly accurate numerical 
bounds we present a method for improving the previously obtained bounds in Section Finally, 
Section is devoted to the derivation of very accurate results for classes of lambda terms which 
have been previously studied in the literature. 
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2 Main results 


In order to count lambda terms of a given size we set up a formal equation which is then translated 
into a functional equation for generating functions. For this we will utilise the symbolic method 
developped in [5]. 

Let us introduce the following atomic classes: the class of zeros Z, the class of successors S, 
the class of abstractions U and the class of applications A. Then the class Coo of lambda terms 
can be described as follows: 


C^ = Seq{S)x Z + U X C^+Ax Cl, (1) 

The number of lambda terms of size n, denoted by Too,n) is |{t G C^o '■ |f| = n}\. Let Loo(z) = 
X]n>ogenerating function associated with C„. Then specification Q gives rise to 
a functional equation for the generating function 

OO 

LooW (2) 

3=0 


Solving 0 we get 


1-z^ 

Cooi^Z^ — 



1-Z>> 


which defines an analytic function in a neighbourhood of z = 0. 


Proposition 1. Let p = RootOf {(1 — z^){l — z'^Y — 4z“+‘^}. Then 


-Loo('^) ^OO ^oo 



+ 0 



(3) 


for some constants a^o > 0, boo < 0 that depend on a, b, c, d. 

Proof. Let f{z) = (1 — z^)(l — — 4z““'''^. Then p is the solution of f{z) = 0. If we compute 

derivative /'(z) = —4(a + — 2cz'^“^(l — z^)(l — z^) — &z^“^(l — z'^)^ we can observe that 

all three terms are negative for 0 < z < 1. Since 0 < p < 1, the function Loo{z) has an algebraic 
singularity of type | which means that its Newton-Puiseux expansion is of the form 

Since Loo{z) is a power series with positive coefficients, we know that Ooo = Loo(p) > 0 and 
feoo < 0. □ 

Corollary 1. The coefficients of Loo{z) satisfy [z”]Loo(z) ~ Cp~'^n~^^‘^, as n ^ oo, where 
C = -boo/{2yCr). 

Let us define the class of m-open lambda terms, denoted Cm, as 


Cm = {t C Coo ■ a prefix of at most m abstractions A is needed to close the term} . 

The number of m-open lambda terms of size n is denoted by Lm,n and the generating function 
associated with the class by Lm{z) = X)n>o ^m.nZ^- Similarly to Coo, the class Cm can be specified, 
and this specification yields the functional equation 

m—1 

Lm{z) = Z^Y. + ^"Lm+l{z) + Z^Lmizf (4) 

3=0 

Note that Lq{z) is the generating function of the set Cq of closed lambda terms. 


3 







Let /Cm = >Coo \ C-m and Km{z) = Loo(-z) — Lm{z). Then using ([^ and Q we obtain 

OO 

Km{z) = z'"^ + Z^'Km+liz) + z’^ Km{z)L^{z) + z‘^ Km{z)Lm{z). 

j^m 


which implies 

Km{z) = 


ya-\-bm 


(1 - z'>)(l - Z‘^{Looiz) + Lm{z))) 1 - Z^{L^{z) + Lm{z)) 


^m+liz'). 


( 5 ) 


( 6 ) 


Note that Km{z) as well as Lm.{z) define analytic functions in a neighbourhood of z = 0. 

Let us state the main theorem of the paper: 

Theorem 1. Let p = RootOf {(1 — z^)(l — Then there exist positive eonstants 

C_ and C (depending on a,b,c,d and m) such that the number of m-open lambda terms of size n 
satisfies 

T iri r^"l T (^\ 

(7) 


liminf&^^lfM > 1 and limsup&t^<L 


n->-oo Cn~2p-'^ 


n->oo Cn 2 p " 


Remark 1. In case of given a,b,c,d and m we can compute numerically such constants C and 
C. This will be done for some of the models mentioned in the introduction. 

Before proving this theorem we will present the key ideas needed for our proof. We introduce 
the class of lambda terms in Cm where the length of each string of successors is bounded by 
a constant integer h. As before, set Lm]n = G C^m ■ |/| = n-| and Lm\z) = ^ 

Then L^m\z) satisfies the functional equation 


T W .,n 

„>Q.l^m,nZ . 


T(h)( + if m</i, 

“ (^) + \im>h. 


( 8 ) 


Notice that for m > h we have a quadratic equation for Lm\z) = L''^\z) that has the solution 


-{h). 


L^h\^) = 


1 - - y'(l - z^)2 - 4z“+‘^^ 

2^^ 


For m < h we have a relation between Lm^ (z) and Li^li(z) which gives rise to a representation of 
{z) in terms of nested radicands (c/. E) after all. Indeed, for m < h we have 


Lt\z) = ^ 



1 - 


(9) 


where 


■l_4^a+dl 


- 2z= 


if m < j < h — 1, 


1 — 2^ 

rj(z) = p - 4z“+‘^i=f4f^ “ 2z^ + 2^^"" ifj = h-l, 

(1 - z^)2 - 4:Z‘^+di^ if j = h. 


One can check that the dominant singularity pm^ of Lm\z) comes from the real root of ^^(z) which 
is closest to the origin and that it is of type 7 {i.e. Lm\z) = ^1 — ^ +0 ( 1 ~ ) 

as z —>■ pm^ for some constants am\b^m depending on m and h). Notice that for all m,k> Owe 
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have pm^ = and that is why we can drop the lower index and write just p^^^ instead of pm'^. 
Moreover, p^^'> > p as well as limft,_>oo P^^'^ = P- 

Let us begin with computing the radii of convergence of the functions Km{z) and Lm{z). For 
the case of binary lambda calculus Lemmas and were already proven in [5]. To extend those 
results to our more general setting, we will use different techniques. 

Lemma 1. For all m> 0 the radius of convergence of Km{z) equals p (the radius of convergence 
of L^{z)). 

Proof. Inspecting (j^ reveals that the key part is — (z)+l (z)) ■ This is the generating function 

of a sequence of combinatorial structures associated with the generating function z^{Lao(z)+Li(z)). 
One can check that we are not in the supercritical sequence schema case {i.e. a singularity of 
considered fraction does not come from the root of its denominator, see [51 PP- 293]) because 
1 — (Loo(p) + Li{p)) > 0. This follows from 

(Loo(p) + U{p)) < 2p‘^LM = 1 - p" < 1. 

The first inequality holds because Loo(p) > Li(p) for alH > 0 and the second one because p > 0. 
Moreover, the radius of convergence of Li{z) is larger than or equal to the radius of convergence 
of Looiz) because £i C Therefore, for all m > 0 the radius of convergence of Km{z) equals p, 
the radius of convergence of Looiz). □ 

Lemma 2. All the functions Lmiz), m > 0, have the same radius of convergence. 

Proof. Let pm denote the radius of convergence of the function Lm (z). From the definition of the 
function L^iz) it is known that for all to > 0 and for all n we have [z"] Lmiz) < [z”] Lm+iiz) and 
therefore pm > Pm+i- Moreover, from Q we know 

m —1 

Lm+liz) = -z“-= ^ Z'^^ + Z-^Lmiz) - Z<^-‘^Lmiz)\ 

3=0 

Notice that pm < 1 because pm < P*-^^ < 1 for > to. Then due to the fact that z~'^Lmiz) — 
z'^~^Lmiz)'^ has radius of convergence bigger or equal pm, we have pm < Pm+i- D 

Lemma 3. For all m > 0 the radius of convergence of Lmiz) equals p. 

Proof. Take , defined in ([^ . Recall that p^"*), pm and p denote the radii of convergence of 
Liff'\z), Lmiz) and Looiz), respectively. Notice that for all m,n > 0 we have [z"]L^^(z) < 
[z^]Lmiz) < [z^] Looiz). Thus we also have p^™^ > Pm > P- Since p^"*^ is the smallest positive 
solution of the equation (1 — z)° — 4z°+‘^ \-zf> ~ have limm->.oo P*-™^ = P, therefore Lemma 2 
implies pm = P- tl 

In the next sections we will present how to obtain an upper and a lower bound for [z"] Lmiz). 
The idea is to construct auxiliary functions satisfying certain inequalities and to use them to 
construct further ones until we have the desired bound. The procedure follows the flowchart 
depicted in Fig. 


3 Upper bound for [z^] L^i^z) 

Notice that for all integers h and to we have Lm'’ C £m- Moreover, for all m,h > 0 there exists 
Um^ such that [z"] Lm\z) = [z"] Lmiz) if n < Um^ and [z”] Lm\z) < [z”] Lmiz) else. We will use 
those properties of Lm\z) in order to derive a lower bound for the asymptotics of [z”] Kmiz). 
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Figure 1: The diagram illustrates the idea how we obtain the upper and the lower bound (in terms 
of the coefficients) for the function L^iz). Starting point is denoted by a blue node, the finish 
nodes are red. 


Note that ([^ corresponds to an equation of the form K-m = -F {ICm, K-m+i, ^m)- Now, de¬ 
fine the new set )Cm^ := T £oo, • From the construction of T and the properties 

of we know that ICm^ C ICm- Let Km^n = ^ ■ 1*1 = Km\z) = X)n >0 Km'!nZ^. 

Then Km\z) satisfies the functional equation 


^ + z‘^K^y{z)L^{z) + z'^K^J^\z)L^J^\z). 


( 10 ) 


In fact, what we did is that we replaced in the application operation every m-open lambda term 
(corresponding to the subterm z^Km{z)Lm{z) of ®) by an m-open lambda ter m wh ere each string 
of successors has bounded length (corresponding to z'^Km{z)Lm\z)). Solving (101 we get 




OO J 

J{b+c) 


n 


l-z'^ (^Looiz) + l\^\z)^ 


— • Sm,oo ('^) ■ 


( 11 ) 


Lemma 4. Let p, Ooc^oo be as in Proposition 0 and Ci = 1/(1 - p'^iaoo + Lf'\p))) and dt 
baop'^/{l — p'^ifloo + l[^\p)))'^. Then Km\z) admits the expansion 




1-^ 

p 


, as z ^ p, 


( 12 ) 


where 


'-'im 


= 


Sm,h-i{p) + R(h) ifm<h, 

_ , 

(l-p'’)(l-p'’+‘=-P‘^(acx=-l-i/'*'(p))) ’ 

EjZL ELm I Ui=m Cfc + RdW ifm<h, 

booP°+’’*”+‘^ 


(^l_p 6 +c_pd (p)) ) ' 


else. 
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with 


Rjh) = 


- 


pa-\-bh-\-c{h—m) 

^-1 

TT r 

(1 - pb) (l - pb+- - pd (uoo + {p)\ 

\\ 11 

j j i—m 

L ^a-\-bh-\-c{h—m)-\-d 

/h-i 

TT 7 

(1 - pb) (l - pb+c - pd (aoo + (p)] 

)) lU 

'' ('^c \ ^ 


\i=m 1 - p^+‘= - p'^ ( Ooo 

+4"^(p) 


Proof. Let us recall that for all m > h we have Lm\z) = L]^ (z). Therefore we can split the 
infinite sum Sm,oo(z) in (111 into the finite one Sm,h-i(z) and the rest Sh,ooiz). 

Case I: m < h. First, consider the finite sum Sm.h-i{z). As in the proof of Lemma we 
identify the key term, show that we are not in the supercritical case, and expand by means of 
Proposition [2 Eventually, this yields 


n 


tta ^ - z:'^iLooiz) + lI^\z)) 


- Cri 


1-^ 


O 


1-^ 


where Cmj = JlL 


Ci and dm ,3 


1 

= Y^i=m % ni=m Cfc- Since ^1 - does neither depend on m 


nor on j and \_^b has poles only on the unit circle, for all m > 0 we have 

Sm,h—l{z') = Cm “f dj^ 


1 - - 
P 


o 


1-^ 


where Cm = J2'j=L and dm = J2’j=L dm,i- 

Now, let us look on the infinite part of the sum in (0. We will refer to its contributions to 
the first and second coefficient of the Newton-Puiseux expansion (121 as the remainders R (h) and 


R,(h), respectively Since for all to > h we have Lm\z) = the sum can be rewritten as 


-(h), 


' h-l 


Sh,oo{z) = 


y.a-\-bh-\-c{h—m) 


ih-zd (^L^{z) + Lf\z)) / (1 - zb) (l - - zd (l^{z) + L^l:\z 


We already know how to handle the product part of this expression, so let us consider the fraction 
-7- 2 -- —Similarly to before, we have to check that the singularity of this 

function does not come from the root of the denominator but from Lao{z) (it cannot come from 
because it has a bigger radius of convergence than Loo{z)). So, we have to show the inequality 

i-pb+^-p'^ (lm + l 

P‘^{Looip) 


(h) 


> 0. But from Loo(p) > L)^ \p) and 0 < < 1 we obtain 

d^h^\p)) — 2p‘^Loo(p) = 1 — < 1 — and hence the desired inequality indeed 

holds. Now, similarly to the previous case we use the Newton-Puiseux expansion of Loo(-z) at p to 
derive an expansion of the infinite part of the sum in (111 and get the asserted result. 

Case II: m > h. This case is easier, because the finite part of the sum in 0 does not exist 
and the other part can be evaluated to a closed form which can be treated as above. □ 


Using the transfer lemmas of |3] (applied to Km\z)) and [z^‘ 
liminf„_,.oo ([z”] Km{z)) ■ V{-\/2)r?l'^p^/d^^ > 1. 


K^\z) < [z'^]Km{z), we get 
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Corollary 2. The number of m-open lambda terms of size n satisfies 


lim sup — 
n^oo Cn~ip~^ 




< 1 where C = 


-'oo 

r(-^) 


4 Lower bound for [z'^] Lm{z) 

The idea behind obtaining a lower bound for [z”] Lm{z) is similar to the one used for the upper 
bound. First we will find an upper bound for [z"] Km{z) using the function 


^ ^ ^ \ L^iz) - K^\z) if TO < i?, 

™ Vooiz) else. 


(13) 


Notice that for all m,h^H^n > 0 we have [z^]Lm{z) < [z"]Lm’^^. Let Cm’^^ denote the 
class of combinatorial structures associated with Lm’^\z) and define the new set '■= 

T >Coo, • From the construction of T and the above properties of Cm'^^ we 

know that fCm fk ■ Let = ||< G : \t\ = n| and Km'^\z) = X)n >0 Km’n'’z'^. 

Then (z) is given by 

OO 

Kt^Hz) = z“ ^ z^^ + z^KiYih^) + z^Ki^’"Hz)Lo.(z) + z‘‘Kt^Hz)L<;^’^\z). 

j^m 


Solving this equation and using (131 we get 

3 


H-l 


A'S'"!-) = TZ?r E n 


j-m 

H-l 


= t _V 

1 - z'> ^ 


yO{b+c) 


1 _ 

=m -L ^ 
3 


^-^oo ('2^) 




ya-\-bH -\-c{H—m) 


TT -_---- 

(2Loo(z)-i^f^W) 


'H-l 

n 


(1 _ (1 _ ^6+c _ 2 z^L^(z)) ^ li 1 (^2Loo{z) - k]^\z)) ] ■ 


(14) 


Lemma 5. Let p be the radius of convergence of the function L^oiz). Then the generating function 
Km’^\z) admits the following expansion 




1-^ 


(15) 


where 


Zh,H) _ P_ 


H-l 


=E n 


1 


i-p'' ^ 

j^m 


H-l 


=m 1 - (SOoo - 


R 


(h,H) , 


J p'^(2b^-dfA 

E E . E. E n 


\ — pb \ — nd- ( 2 n — ^ — nd (20 — 

j—m i—m ^ • \ j i—m ^ ^ \ j 


+ Rjh,H), 



















Oooi&oo and come from the expansion of Loo{z) andK[^\z), respectively, at p (see Propo¬ 

sition^ and the proof of Lemma^ and 


R {h.H) = 


R.(h,H) = 


p' 


fL-\-hH. -\-c{H—m) 


' H-1 

n 


(1 - pV (1 - - 2p^a^) 1 _ ( 2 „^ _ J ■ 

.a-\-bH-\-c(H—m)-\-d f H—1 ^ 

n 


(1 - p^} (1 - - 2p‘‘a^) 11 1 _ (■ 2 „^ _ 

X 



As in the previous section, we apply the the transfer lemmas of a to Km’^\z) and use 
[z"] Km’^\z) > [z”] Km{z) to arrive at the following result: 


Corollary 3. The number of m-open lambda terms of size n satisfies 


liminf > 1 where C = ^21^ 


n->-oo Qji 2 p-ii 


r(-i) 


5 Improvement of the bounds 

The bounding functions L^^z) = Lao{z) — Km’^\z) and L^^z) = L,x,{z) — Km\z) derived in the 
previous sections can be used in a straightforward way to compute numerical values for C and C, 
given concrete values for a, b, c, d. If we choose h and H big enough, then this will give us proper 
bounds. But in practice they still leave a large gap, at least for values h and H that allow us to 
perform the computation within a few hours on a standard PC. For instance, in case of the natural 
counting for ft, = iJ = 15 we get « 0.00404525... and (7*-"'“*^ ~ 0.18086721.... 

We show a simple way to improve C and C. Let us introduce the functions L^^j,j-{z) and 
L^’^^(z), defined by the following equations: 


L 




m+l,M V 


W + if m<M, 


if TO = M, 


^m,M 


(z) = 


+ z‘^l'^^]’^{zY if to < M, 

- K^m'^\z) 


if TO = M. 

These two function admit a representation in terms of nested radicals which is similar to ([^ : 


-W 


,,M 


(z) = 




2zd 

_ 

^m,M — 

1 f 


l-4z‘^+<‘ -2z‘=+2z<= 


y...yi- 4 z°+d -2z"+2z‘=A/roo(z)-J^l?^z)j 


1 - 


l-d^a+d -2z^-\-2z^ 




So, in and {z) we are using exact expressions for Lm{z) up to some constant M 

and then replace Lm{z) by a function that is its upper and lower bound, respectively. Numerical 
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results for some previously studied notions of size (see Sections 6.1 and 6.21 reveal a significant 
improvement in closing the gap between the constants (7, C obtained by utilizing the functions 


-(h) 


„M 


(2:)- 


6 Results for some previously introduced notions of size 

6.1 Natural counting 

Lemma 6. The following bounds hold 


lim inf 

n—>-oo 




> 1 


and 


lim sup 

n—>-oo 


[z"] 

y^(nat) _3 _ 

C n ‘2 p 


< 1 


(16) 


where p = RootOfj—1 + 3x + « 0.295598... and ^("“0 are computable constants 

with numerical values « 0.00404525 ... and (7*'"°*^ « 0.18086721.... 


Proof. For ’natural counting’ we have a = b = c = d = l. From Theorem we know that the 

radius of convergence of equals p = RootOf { —1 + 3a: + + x^} « 0.295598 We can 

also easily get Loo(z) ^ a^a + 6oo\/l — z/p with 


Q-oo 


P 

2p 


1.19149 ... and boo 


1 /1 + p + 

;^Y 2^ 


2.15093..., 


and from the transfer lemmas of m we obtain [z”] Loo (- 2 ) ~ "/r(—1/2) ~ (0.606767.. .)• 

n“5 (3.38298 ...)", as n —>■ oo. 


h,H 

^0 


^0 

AhM) 

1 

0.855448 

-1.153959 

1.086200 

-3.803686 

2 

0.898032 

-1.313246 

0.979519 

-2.581823 

3 

0.917305 

-1.397536 

0.958215 

-2.324953 

4 

0.927248 

-1.444672 

0.950295 

-2.236290 

5 

0.932849 

-1.472308 

0.946185 

-2.192353 

6 

0.936128 

-1.488826 

0.943824 

-2.167379 

7 

0.938055 

-1.498647 

0.942443 

-2.152790 

8 

0.939174 

-1.504385 

0.941643 

-2.144335 

9 

0.939813 

-1.507673 

0.941187 

-2.139511 

10 

0.940172 

-1.509525 

0.940931 

-2.136799 

11 

0.940372 

-1.510556 

0.940788 

-2.135291 

12 

0.940482 

-1.511125 

0.940710 

-2.134460 

13 

0.940543 

-1.511438 

0.940667 

-2.134004 

14 

0.940576 

-1.511608 

0.940643 

-2.133755 

15 

0.940594 

-1.511701 

0.940630 

-2.133619 


Table 1: Numbers rounded up to 7 digits. 


Since we are most interested in the enumeration of closed lambda terms, we examine the mul¬ 
tiplicative constants in the leading term of the asymptotical lower and upper bound for [z"] Lq(z). 
From the formulas in Lemmas andwe have computed the values for Cq^\ d^'^ and Cq^’^\ d^’^^ 
for different constants h and H (see Table . 
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0.18086721 0.00404525 
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n‘2 p-" 



Figure 2: Plot of the exact value of 
15 and 10 < n < 150. 


_ 3 
n 2 


and computed lower and upper bound for h = H = 


As expected, the bigger h and H are, the more accnrate is the bound we get (in this case 
for [z'^] Ko{z)). Taking the values of and for h = H = 15, Corollaries and yield 

Q(nat) ^ 0.00404525... and ~ 0.18086721.... Notice that nsing the valnes of to 

compute gives non-trivial values only for > 7 (for 1 < h < 7 we get negative numbers 


because in this case we have 




> l&ool). 


□ 


Fignre|^illnstrates the bonnds we obtained and the exact values of the coefficients [z"] (z) 

for 10 < n < 150. Applying the approach discnssed in Section for h = H = M = 13 we get the 
following improvement. 


Lemma 7. The following bounds hold 


lim inf 

n—>-oo 


[z"]4”“‘^(z) 


> 1 and lim sup 

n—>-oo 


[z"] 4"“‘^(z) 


— (nat) 

C n 


3 

2/9 


— n 


< 1 


(17) 


/ — (nat) 

where p = RootOfj —1 + 3a; + ~ 0.295598 ... and Q}'^^ \ C are computable constants 

, ={nat) 

with numerical values « 0.07790995266 ... and C « 0.07790998229 .... 

Figurej^illustrates how the improvement discnssed in Section [fallows to reduce the gap between 

, . —(nat) 

the constants ',C for the lower and the upper bound. 


6.2 Binary lambda calculus 

Lemma 8. The following bounds hold 


lim inf 

n—>-oo 




> 1 


and 


lim sup 

n—>-oo 


[z’ 


L<“”'(=) 


— (bin) 

C n 


2p " 


< 1 


(18) 
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c 



. — (bin) 

where p = RootOf{—1 + a; + 2x^ — 2x^ + 3a:^ + x^} « 0.509308 ... and C are computable 

, , =(bin) 

constants with numerical values C} « 0.01252417... and C « 0.01254593 .... 

In order to proof this lemma it is enough to recall that in case of binary lambda calculus the size 
defining constants are 6=1 and a = c = d = 2. Then we used the functions Lq\^^{z), to 

obtain the numerical constants stated in Lemma [HI 
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Appendix 

Below we present the proof of Lemma 


Proof. Throughout the proof we will make use of Equation (141. Let us first focus on the finite sum: 

^-L-The key part of it - ^^is the 


1—z^ 




- 7 - Y7 -.—The key part of it -?- 77 -,—v 

l-z‘^(2L^(z)-Kl^\z)] l-z<ihL^(z)-Kl’'fzU 

generating function of a sequence of structures enumerated by the denominator z'^ {^Loo (z) — (z) 

Like in case of Equation (11), one can check that we are not in the supercritical sequence schema 
case because 

l-p'" {2LM - k\'^\p)) >0 

{2LM - K^t\p)) < 2p"Loo(p) = 1 - < 1. 

Second inequality holds because k\^\p) > 0 for all i > 0. The last inequality holds because we 
have p > 0 for all a, b, c, d that satisfies Assumption]^ Therefore, for all i from the Newton-Puiseux 
expansion of Lao{z) and K^^\z) at singularity p (see Proposition and Lemma we have : 


1 


1-z‘i (^2L^{z) - k['^\z) 


= Ci 


d,; 1 - 


+ 01 - 




where c,; = -— - - 71 :^, di = ^^- ' L ^ 2 . Moreover, we have 


(19) 


fr_i_ 

lX\-z<^{:iL^{z)-Kf\z)) 


^m,j “f 


1 - 


+ 01 - 


where c^j = ^Ili=m^k- Since (^1 - " does not depend on m and j 

and \_^b have poles only on the unit circle, we have for all m > 0 the following 

a—cm 3 

Y^b E n 


1 


(2L^{z) - K^,^\z) 
where c" — „j(6+c)a . d' _ 

wneie Cm — i_pb 2^j=mP Cmj,Um— i_pb Umj- 


1 - - 
p 


0 ( 1 -- 
p 


In case of the remainder: 


a + bH + c(H-m.) 


(l-z‘^)(l-z^+=-2z‘*Loo(z)) 


H-1 


nr= 


- 7 - - - 777 —V I we proceed 

l-z‘i( 2 Loo(z)-if<'‘>(z)) J 

similarly to the proof of Lemma We already know how to handle the product part of this 
expression, so let us consider the fraction: T-z-^f(i-z‘’+‘=- 2 z‘^L —(!))■ Similarly to before, we have 
to check that the singularity of this function does not come from the root of the denominator but 
from the function Lco{z). So the following inequality has to hold: 

l-p''+^-2p^L^(p)>0 
2p‘^Loo(p) = 1 - p" < 1 - 

Second inequality holds because we have 0 < p^^'^ < p'^ for all o, 6 , c, d that satisfies Assumption]^ 
Now, similarly to the previous case using the Newton-Puiseux expansion of Laoiz) at p we can 
derive expansion of this part of Equation (14). □ 
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